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METHOD AND SYSTEM TO INVESTIGATE A 
COMPLEX CHEMICAL SPACE 



This invention was made with government support under Contract No. 
70NANB9H3038 awarded by NIST. The government may have certain rights to the 
invention. 

BACKGROUND OF THE INVENTION 

[0001] The present invention relates to a combinatorial high throughput 

screening method and system. 

[0002] Combinatorial organic synthesis (COS) is a high throughput screening 

method that was developed for pharmaceuticals. COS uses systematic and repetitive 
synthesis to produce diverse molecular entities formed from sets of chemical "building 
blocks." As with traditional research, COS relies on experimental synthesis 
methodology. However instead of synthesizing a single compound, COS exploits 
automation and miniaturization to produce large libraries of compounds through 
successive stages, each of which produces a chemical modification of an existing 
molecule of a preceding stage. Libraries are physical, trackable collections of samples 
resulting from a definable set of the COS process or reaction steps. The libraries 
comprise compounds that can be screened for various activities. 

[0003] Combinatorial high throughput screening (CHTS) is an HTS method 

that incorporates characteristics of COS. The CHTS methodology is marked by the 
search for high order synergies and effects of complex combinations of experimental 
variables through the use of large arrays in which multiple factors can be varied 
through multiple levels. Factors of an experiment can be varied within an array 
(typically formulation variables) and between an array and a condition (both 
formulation and processing variables). Results from the CHTS experiment can be 
used to compare properties of the products in order to discover "leads" - formulations 
and/or processing conditions that indicate commercial potential. 
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[0004] The steps of a CHTS methodology can be broken down into generic 

operations including selecting chemicals to be used in an experiment; introducing the 
chemicals into a formulation system (typically by weighing and dissolving to form 
stock solutions), combining aliquots of the solutions into formulations or mixtures in 
a geometrical array (typically by the use of a pipetting robot); processing the array of 
chemical combinations into products; and evaluating the products to produce results. 

[0005] Typically, CHTS methodology is characterized by parallel reactions at 

a micro scale. In one aspect, CHTS can be described as a method comprising (A) an 
iteration of steps of (i) selecting a set of reactants; (ii) reacting the set and (iii) 
evaluating a set of products of the reacting step and (B) repeating the iteration of steps 
(i), (ii) and (iii) wherein a successive set of reactants selected for a step (i) is chosen as 
a result of an evaluating step (iii) of a preceding iteration. 

[0006] Results from an experiment can be evaluated by aid of mathematical 

models as taught by G.E.P. Box and N. R. Draper, Empirical Model-Building and 
Response Surfaces, John Wiley and Sons, NY, 1987, p 20-22. The models can be 
used to find an approximation, typically a polynomial, to an unknown underlying 
theoretical function. For example, Taylor's Series expansion is a polynomial that can 
provide a valuable approximation of first or second order experimental interactions. 

[0007] However, the study of catalyzed chemical reactions by CHTS involves 

the investigation of a complex experimental space characterized by multiple 
qualitative and quantitative factor levels. Typically, the interactions of a catalyzed 
chemical reaction such as a carbonylation reaction can involve interactions of an order 
of 6 or 9 or greater. While Taylor expansion approximation can be effectively applied 
to analyze first or second order interactions, it is useless to study CHTS results from a 
complex catalyzed chemical reaction. 

[0008] Another problem is that catalyzed chemical reactions can be 

unpredictable. Well-known protocols in one area of chemistry cannot be applied to 
another area with assurance of success. For example, U.S. Pat. 6143914 shows that 
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some combinations of various metals unexpectedly increase a carbonylation catalyst 
turnover number (TON) and other related combinations do not. "Due to the 
complicated mechanistic nature of many transition metal based catalysts, structure - 
activity relationships are often unpredictable, leaving empirical exploration and 
serendipity the most common routes to discovery." J. Tian & G. W. Coates, Angew. 
Chem Int. Ed. 2000, 39, p 3626. This high degree of irregularity and unpredictability 
is illustrated in FIG.l. 

[0009] There is a need for a methodology to examine the complex higher 

n order and unpredictable interactions of a CHTS catalyzed chemical reaction 
*B experiment that cannot be examined by Taylor Series expansion or other standard 

iy methodology. 



BRIEF SUMMARY OF THE INVENTION 



[0010] The invention provides a particularly well-suited experimental 

methodology to investigate multiple and complex interactions of a catalyzed chemical 
fi reaction that involves both qualitative and quantitative factor levels. According to the 
S3 invention, an experimental space of a catalyzed chemical reaction is defined to 

represent at least three factor interactions, a CHTS method is effected on the catalyzed 
chemical experimental space to produce results and results are analyzed according to 
matrix algebra to select a best case set of factor levels from the catalyzed experimental 
space. 

[00 1 1 ] In another embodiment, a CHTS experiment is conducted on a 

complex experimental space comprising qualitative and quantitative factors to 
produce first data results, (B) the first data results are analyzed according to matrix 
algebra, (C) a standard deviation of the analyzed results is defined, (D) data results 
that positively exceed the standard deviation are selected, (E) a next experimental 
space is defined according to the selected data results and (F) steps (A) through (E) 
are reiterated on the next experimental space until data results selected in step (D) 
represent satisfactory leads. 
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[0012] In yet another embodiment, a system for investigating a catalyzed 

experimental space, comprises a reactor for effecting a CHTS method on the catalyzed 
chemical experimental space to produce results and a programmed controller to 
analyze the results according to matrix algebra to select a best case set of factor levels 
from the catalyzed experimental space. 

BRIEF DESCRIPTION OF THE DRAWINGS 



[0013] FIG. 1 is a graphic representation of a complex catalyzed chemical 

experimental space; 

[0014] FIG. 2 is a schematic representation of a CHTS system; 

[0015] FIG. 3 shows representations of vector matrix forms; and 

[00'6] FIG. 4 is a normal probability plot. 

DETAILED DESCRIPTION OF THE INVENTION 

[0017] A purpose of a CHTS experiment is to locate anomalies that represent 



high-value "leads." Leads are results that identify candidate factors and levels for a 
commercial process. According to the invention, matrix algebra is combined with 
CHTS to examine the very high-level interactions of a chemical catalytic reaction. In 
the CHTS/matrix algebra approach, a CHTS experiment in many factors can be run 
and then modeled to generate all main, 2 nd , 3 rd , 4 th and even higher interactions. 
Suitably, the results can be represented according to a general linear model (GLM). A 
statistical analysis of the represented results determines high-order interaction leads. 

[0018] In a typical CHTS method, a multiplicity of tagged reactants is 

subjected to an iteration of steps of (A) (i) simultaneously reacting the reactants, (ii) 
identifying a multiplicity of tagged products of the reaction and (B) evaluating the 
identified products after completion of a single or repeated iteration (A). A typical 
CHTS can utilize advanced automated, robotic, computerized and controlled loading, 
reacting and evaluating procedures. 



RD-28250 




[0019] The results from the CHTS experiment are analyzed using matrix 

algebra to extract combinations of the experiment interactions. A mathematical 
matrix is a representation of real numbers in a rectangular array. Matrices are 
important tools for expressing and discussing problems that may involve complex 
data sets. Matrix analysis is a multivariant methodology for expressing and 
manipulating these kinds of data and for solving problems posed by the data. Matrix 
operations can include representation (posing or modeling data in a matrix 
representation), addition, subtraction, scalar multiplication, matrix multiplication, 
multiplication by inverse, transposition (interchanging rows and columns) and 
distribution (assigning a probability value). Matrices can be manipulated to produce a 
sum, difference, scalar multiple, matrix multiple, product or transpose. 

[0020] In the present invention, the matrices are representations of CHTS 

results in a rectangular array. The runs from the CHTS experiment provide sets of 
results or v's, one for each run, each correlated with a set of levels of factors, x t . Each 
run y is associated with an error e. Each of the factors or interactions is associated 
with a coefficient p. These elements (x's, their interactions and /s) can be 
represented in vector/matrix form as shown in FIG. 3, wherein levels of factors and 
interactions form a rectangular array or matrix (20) of scalar values X. Further in FIG. 
3, y's, (?s and <?'s are represented in single column matrices (10), (30), and (40). A 
matrix estimation equation of the system can be as follows: 

y = Xp + e (I) 

[0021] where X is a matrix of factor and interaction levels in the experiment, y 

is a matrix of experimental results, p is effects (both main effects and interaction 
effects), every e, in e has the same variance a 2 , error terms e, arise from a normal 
distribution and expected value E(y) (definition: the value of the response if no error 
is present) of the response is E(y) = Xp. 

[0022] The method involves solving the above matrix estimation equation (I), 

according to the relationship: 
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P = (X'X)' 1 X , y (ID 

[0023] where superscript ' indicates a transpose of a matrix (in which each 

row becomes a column and each column a row). The superscript _1 indicates an 
inverse function of a matrix. Thus for any square matrix A, a relation can be defined 
as AA" 1 = I, where I is the identity matrix 50 shown in the FIG. 3 model. 

[0024] Accordingly, results can be assembled as an n x 1 vector y and factor 

level values can be assembled into annxfc+1 matrix X with l's as designations in a 
first column and each other column containing the coded factor level values ( +l's and 
-l's representing the extents of the values of the factors and interactions). Matrix 
equation (I) is then solved for effects parameters p. 

[0025] The effects parameters of matrix p are then examined for statistical 

significance. The null hypothesis can be applied in this examination. The null 
hypothesis is that all of the effects observed in the experiment are caused simply by 
random processes. If this is correct, the effects will fit to a normal distribution and 
form a relatively straight line in a probability plot. In FIG. 4; E100 is a straight line 
representing a results approximation. The line is flanked by dashed lines denoting 
multiples of standard deviations. A desired standard deviation can be selected by an 
experimenter for the experiment. Any effects that fall off the line by more than the 
standard deviation can be interpreted to have been caused by nonrandom processes, as 
taught by D. Montgomery, Design and Analysis of Experiments, 3 rd Ed., John Wiley, 
1991, NY, p 99. Effects that positively exceed the deviation can represent 
combinations that are failures or combinations that provide synergistic improvement, 
i.e., leads. 

[0026] In a preferred embodiment, results from the CHTS method are 

analyzed by matrix algebra by steps of (A) representing the results as an n x 1 matrix y 
where n = a number of factor level combinations in the experiment, (B) representing 
extents of the factor level combinations in an n x n matrix X, (C) solving n 
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simultaneous equations represented by the matrices according to matrix algebra to 
form a results matrix P and (D) examining the results matrix (3 to identify effects 
outside a standard deviation. 

[0027] The step (B) can comprise coding extents of the factor level 

combinations as a +1 or -1 and representing the coded extents as the n x 1 matrix y. 
The step (C) can comprise (i) transposing matrix X to form matrix X', (ii) 
postmultiplying X' by X to generate a matrix and (iii) postmultiplying the generated 
matrix by y to form the results matrix (3. The step (D) can comprise (i) representing 
the results matrix p as a normal probability plot, defining a standard deviation for 
results of the plot and (iii) identifying positive interactions outside of the standard 
deviation. The standard deviation can represent a probability that a result deviation 
from the standard is random and that positive interactions can be identified outside of 
the deviation. In one embodiment, the probability can be established at 95 percent or 
better to define an experimental space for a commercial process or the probability can 
be established at 99.7 percent or better to define a best set of factor levels as leads for 
a commercial process. 

[0028] I n one embodiment, the invention is applied to screen for a catalyst to 

prepare a diaryl carbonate by carbonylation. Diaryl carbonates such as diphenyl 
carbonate can be prepared by reaction of hydroxyaromatic compounds such as phenol 
with oxygen and carbon monoxide in the presence of a catalyst composition 
comprising a Group VIHB metal such as palladium or a compound thereof, a bromide 
source such as a quaternary ammonium or hexaalkylguanidinium bromide and a 
polyaniline in partially oxidized and partially reduced form. 

[0029] Various methods for the preparation of diaryl carbonates by a 

carbonylation reaction of hydroxyaromatic compounds with carbon monoxide and 
oxygen have been disclosed. The carbonylation reaction requires a rather complex 
catalyst. Reference is made, for example, to Chaudhari et al., U.S. Pat. 5,917,077. 
The catalyst compositions described therein comprise a Group VlllB metal (i.e., a 
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metal selected from the group consisting of ruthenium, rhodium, palladium, osmium, 
iridium and platinum) or a complex thereof. 

[0030] The catalyst material also includes a bromide source. This may be a 

quaternary ammonium or quaternary phosphonium bromide or a 
hexaalkylguanidinium bromide. The guanidinium salts are often preferred; they 
include the V, T-bis(pentaalkylguanidinium)alkane salts. Salts in which the alkyl 
groups contain 2-6 carbon atoms and especially tetra-n-butylammonium bromide and 
hexaethylguanidinium bromide are particularly preferred. 

[0031] Other catalytic constituents are necessary in accordance with 

Chaudhari et al. The constituents include inorganic cocatalysts, typically complexes 
of cobalt(II) salts with organic compounds capable of forming complexes, especially 
pentadentate complexes. Illustrative organic compounds of this type are nitrogen- 
heterocyclic compounds including pyridines, bipyridines, terpyridines, quinolines, 
isoquinolines and biquinolines; aliphatic polyamines such as ethylenediamine and 
tetraalkylethylenediamines; crown ethers; aromatic or aliphatic amine ethers such as 
cryptanes; and Schiff bases. The especially preferred inorganic cocatalyst in many 
instances is a cobalt(H) complex with bis-3-(salicylalamino)propylmethylamine. 

[0032] Organic cocatalysts may be present. These cocatalysts include various 

terpyridine, phenanthroline, quinoline and isoquinoline compounds including 
2,2':6 , ,2"-terpyridine, 4-methylthio-2,2':6',2 H -terpyridine and 2,2':6',2"-terpyridine N- 
oxide, 1 , 10-phenanthroline, 2,4,7,8-tetramethyl- 1 , 10-phenanthroline, 4,7-diphenyl- 
1,10, phenanthroline and 3,4,7,8-tetramethy-l, 10-phenanthroline. The terpyridines 
and especially 2,2':6',2 M -terpyridine are preferred. 

[0033] Another catalyst constituent is a polyaniline in partially oxidized and 

partially reduced form. 

[0034] Any hydroxyaromatic compound may be employed. 

Monohydroxyaromatic compounds, such as phenol, the cresols, the xylenols and p- 
cumylphenol are preferred with phenol being most preferred. The method may be 
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employed with dihydroxyaromatic compounds such as resorcinol, hydroquinone and 
2,2-bis(4-hydroxyphenyl)propane or "bisphenol A," whereupon the products are 
polycarbonates. 

[0035] Other reagents in the carbonylation process are oxygen and carbon 

monoxide, which react with the phenol to form the desired diaryl carbonate. 

[0036] These and other features will become apparent from FIG. 2 and the 

following detailed discussion, which by way of example without limitation describe 
preferred embodiments of the present invention. 

[0037] FIG. 2 is a schematic representation of a system 10 for CHTS 

according to the invention. FIG. 2 shows system 10 including dispensing assembly 
12, reactor 14, detector 16 and controller 18. Further shown, is X-Y-Z robotic 
positioning stage 20, which supports array plate 22 with wells 24. The dispensing 
assembly 12 includes a battery of pipettes 26 that are controlled by controller 18. X- 
Y-Z robotic positioning stage 20 is controlled by controller 18 to position wells 24 of 
the array plate 22 beneath displacement pipettes 26 for delivery of test solutions from 
reservoirs 28. 

[0038] Controller 18 controls aspiration of precursor solution into the battery 

of pipettes 26 and sequential positioning of the wells 24 of array plate 22 so that a 
prescribed stoichiometry and/or composition of reactant and/or catalyst can be 
delivered to the wells 24. By coordinating activation of the pipettes 26 and movement 
of plate 22 on the robotic X-Y-Z stage 20, a library of materials can be generated in a 
two-dimensional array for use in the CHTS method. Also, the controller 18 can be 
used to control sequence of charging of sample to reactor 14 and to control operation 
of the reactor 14 and the detector 16. Controller 18 can be a computer, processor, 
microprocessor or the like. 

[0039] An experimental space definition defines the contents of the wells 24 

for the CHTS method. The space can be defined according to any design that results 
in a representation of at least three factor interactions. Suitable designs include 
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fractional factorial design, Latin square design, Plackett-Burman design or Taguchi 
design. Preferably the design results in a representation of all interactions and 
preferably, the design is an orthogonal design such as a full factorial design. The 
design can be embodied as an algorithm or program resident in controller 18. 

[0040] Controller 18 controls the sequence of charging array plate 22 into the 

reactor 14, which is synchronized with operation of detector 16. Detector 16 detects 
products of reaction in the wells 24 of an array plate 22 after reaction in reactor 14. 
Detector 16 can utilize chromotography, infra red spectroscopy, mass spectroscopy, 
laser mass spectroscopy, microspectroscopy, NMR or the like to determine the 
constituency of each reaction product. The controller 18 uses data on the sample 
charged by the pipettes 26 and on the constituency of reaction product for each sample 
from detector 16 to correlate a detected product with at least one varying parameter of 
reaction. Additionally, an algorithm or program can be resident in the controller 18 to 
represent the CHTS results according to a matrix form and to analyze the represented 
results by matrix algebra to determine leads. 

[0041] As an example, if the method and system of FIG. 1 is applied to study a 

carbonylation catalyst and/or to determine optimum carbonylation reaction conditions, 
the detector 16 analyzes the contents of the well for carbonylated product. In this 
case, the detector 16 can use Raman spectroscopy. The Raman peak is integrated 
using the analyzer electronics and the resulting data can be stored in the controller 18. 
Other analytical methods may be used - for example, Infrared spectrometry, mass 
spectrometry, headspace gas-liquid chromatography and fluorescence detection. 

[0042] A method of screening complex catalyzed chemical reactions can be 

conducted in the FIG. 2 system 10. According to the method, catalyzed formulations 
are prepared according to any suitable procedure. For example, one procedure 
produces a homogeneous chemical reaction utilizing multiphase reactant systems. In 
this procedure, a formulation is prepared that represents a first reactant system that is 
at least partially embodied in a liquid. Each formulation is loaded as a thin film to a 
respective well 24 of the array plate 22 and the plate 22 is charged into reactor 14. 
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During the subsequent reaction, the liquid of the first reactant system embodied is 
contacted with a second reactant system at least partially embodied in a gas. The 
liquid forms a film having a thickness sufficient to allow the reaction rate of the 
reaction to be essentially independent of the mass transfer rate of the second reactant 
system into the liquid. 

[0043] The method herein described can be used with any suitable catalyzed 

chemical reactant system. For example, the system and method herein can be used for 
determining a method for producing diphenyl carbonate (DPC). Diphenyl carbonate 
(DPC) is useful, inter alia, as an intermediate in the preparation of polycarbonates. 
One method for producing DPC involves the carbonylation of a hydroxyaromatic 
compound (e.g., phenol) in the presence of a catalyst system. A carbonylation catalyst 
system typically includes a Group Vffl B metal (e.g., palladium), a halide composition 
and a combination of inorganic co-catalysts (IOCCs). 

[0044] Generally, testing of new catalyst systems has been accomplished at 

macro-scale and, because the mechanism of this carbonylation reaction is not fully 
understood, the identity of additional effective IOCCs has eluded practitioners. An 
embodiment of the present invention allows a homogeneous carbonylation reaction to 
be carried out in parallel with various potential catalyst systems and, consequently, 
this embodiment can be used to identify effective IOCCs for the carbonylation of 
phenol. 

[0045] The following Example is illustrative and should not be construed as a 

limitation on the scope of the claims unless a limitation is specifically recited. 

EXAMPLE 

[0046] This EXAMPLE illustrates an identification of an active and selective 

catalyst for the production of aromatic carbonates. The procedure identifies the best 
catalyst from a complex chemical space, where the chemical space is defined as an 
assemblage of all possible experimental conditions defined by a set of variable 
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parameters such as formulation ingredient identity or amount or process parameter 
such as reaction time, temperature, or pressure. 

[0047] The chemical space consists of the following TABLE 1 chemical factor 

levels and TABLE 2 processing factor levels: 



Factor 
Primary Catalyst 



Metal Cocatalyst 



TABLE 1 

Level Le vel 
Ru(acac)3 All at 25 ppm 

Pt(acac)2 

Mn(acac)2 150 and 1500 ppm 

Fe(acac)3 



Cosolvent Dimethylformamide All at 10% 

(DMFA), Tetrahydrofuran 
(THF) 

Anion Cocatalyst CT, Br. (as All at 5000 ppm 

hexamethylguanadinium 
salts) 

TABLE 2 

Factor Level . _ , _ . ... 

Pressure 1000 psi, 1500 psi (8% Oxygen m Carbon Monoxide) 

Temperature 100 C, 120 C 
[0048] The system has seven factors, each at two levels. There are 2 7 = 128 

possible combinations of these levels. The experiment is set up according to a full 
factorial design with 128 runs as shown in TABLE 3. In the experiment, catalyzed 
mixtures are made up in phenol solvent using the concentrations of each component 
as given in the rows of TABLE 3. The total volume of each catalyzed mixture is 1.0 
ml. From each mixture, a 25 microliter aliquot is dispensed into a 2 ml reaction vial, 
forming a film on the bottom. The vials are grouped in array plates by process 
conditions (as specified in the Pressure and Temperature columns in the table) and 
each array plate is loaded into a high pressure autoclave and subjected to the reaction 
conditions specified. At the end of the reaction time, the reactor is cooled and 
depressurized and the contents of each vial are analyzed for diphenyl carbonate 
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product using a gas chromatographic method. A turnover number (TON) for each 
reaction is calculated as mols of diphenylcarbonate/mols of primary catalyst. The 
results are given in the TON column of TABLE 3. 
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TABLE 3 



C: 

A: B: Metal 

Primary Metal Cocataly 

Catalyst Cocatalyst Amount 

Pt Mn 150 

Ru Mn 150 

Pt Fe 150 

Ru Fe 150 

Pt Mn 1500 

Ru Mn 1500 

Pt Fe 1500 

Q Ru Fe 1500 

»3 pt Mn 150 

"'3 R U Mn 150 

W p t Fe 150 

^ Ru Fe 150 

£ Pt Mn 1500 

H Ru Mn 1500 

T Pt Fe 1500 

p Ru Fe 1500 

fi Pt Mn 150 

1U Ru Mn 150 

Si Pt Fe 150 

Q Ru Fe 150 

\* Pt Mn 1500 

Ru Mn 1500 

Pt Fe 1500 

Ru Fe 1500 

Pt Mn 150 

Ru Mn 150 

Pt Fe 150 

Ru Fe 150 

Pt Mn 1500 

Ru Mn 1500 

Pt Fe 1500 

Ru Fe 1500 
Pt Mn 150 

Ru Mn 150 
Pt Fe 150 

Ru Fe 150 
Pt Mn 1500 
Ru Mn 1500 
Pt Fe 1500 



D: 


Anion 


F: 


G: 




Cosolvent 


Cocatalyst Pressure 


Temp. 


TON 


DMFA 


Br 


1000 


100 


3340 


DMFA 


Br 


1000 


100 


3470 


DMFA 


Br 


1000 


100 


2360 


DMFA 


Br 


1000 


100 


2260 


DMFA 


Br 


1000 


100 


2310 


DMFA 


Br 


1000 


100 


2060 


DMFA 


Br 


1000 


100 


3030 


DMFA 


Br 


1000 


100 


3200 


THF 


Br 


1000 


100 


2430 


THF 


Br 


1000 


100 


2270 


THF 


Br 


1000 


100 


2910 


THF 


Br 


1000 


100 


3160 


THF 


Br 


1000 


100 


3270 


THF 


Br 


1000 


100 


3030 


THF 


Br 


1000 


100 


2260 


THF 


Br 


1000 


100 


2470 


DMFA 


CI 


1000 


100 


3040 


DMFA 


CI 


1000 


100 


3340 


DMFA 


CI 


1000 


100 


2030 


DMFA 


CI 


1000 


100 


1860 


DMFA 


CI 


1000 


100 


2200 


DMFA 


CI 


1000 


100 


1920 


DMFA 


CI 


1000 


100 


3290 


DMFA 


CI 


1000 


100 


2910 


THF 


CI 


1000 


100 


2260 


THF 


CI 


1000 


100 


2410 


THF 


CI 


1000 


100 


3260 


THF 


CI 


1000 


100 


3200 


THF 


CI 


1000 


100 


3360 


THF 


CI 


1000 


100 


3090 


THF 


CI 


1000 


100 


2320 


THF 


CI 


1000 


100 


2320 


DMFA 


Br 


1500 


100 


3230 


DMFA 


Br 


1500 


100 


3710 


DMFA 


Br 


1500 


100 


2140 


DMFA 


Br 


1500 


100 


2500 


DMFA 


Br 


1500 


100 


2490 


DMFA 


Br 


1500 


100 


2230 


DMFA 


Br 


1500 


100 


3070 
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Ru 


Fe 


1500 


DMFA 


Br 


Pt 


Mn 


150 


THF 


Br 


Ru 


Mn 


150 


THF 


Br 


Pt 


Fe 


150 


THF 


Br 


Ru 


Fe 


150 


THF 


Br 


Pt 


Mn 


1500 


THF 


Br 


Ru 


Mn 


1500 


THF 


Br 


Pt 


Fe 


1500 


THF 


Br 


Ru 


Fe 


1500 


THF 


Br 


Pt 


Mn 


150 


DMFA 


CI 


Ru 


Mn 


150 


DMFA 


CI 


Pt 


Fe 


150 


DMFA 


CI 


Ru 


Fe 


150 


DMFA 


Cl 


Pt 


Mn 


1500 


DMFA 


CI 


Ru 


Mn 


1500 


DMFA 


CI 


Pt 


Fe 


1500 


DMFA 


CI 


Ru 


Fe 


1500 


DMFA 


CI 


Pt 


Mn 


150 


THF 


CI 


Ru 


Mn 


150 


THF 


CI 


Pt 


Fe 


150 


THF 


CI 


Ru 


Fe 


150 


THF 


CI 


Pt 


Mn 


1500 


THF 


CI 


Ru 


Mn 


1500 


THF 


Cl 


Pt 


Fe 


1500 


THF 


CI 


Ru 


Fe 


1500 


THF 


Cl 


Pt 


Mn 


150 


DMFA 


Br 


Ru 


Mn 


150 


DMFA 


Br 


Pt 


Fe 


150 


DMFA 


Br 


Ru 


Fe 


150 


DMFA 


Br 


Pt 


Mn 


1500 


DMFA 


Br 


Ru 


Mn 


1500 


DMFA 


Br 


Pt 


Fe 


1500 


DMFA 


Br 


Ru 


Fe 


1500 


DMFA 


Br 


Pt 


Mn 


150 


THF 


Br 


Ru 


Mn 


150 


THF 


Br 


Pt 


Fe 


150 


THF 


Br 


Ru 


Fe 


150 


THF 


Br 


Pt 


Mn 


1500 


THF 


JbJr 


Ru 


Mn 


1500 


THF 


Br 


Pt 


Fe 


1500 


THF 


Br 


Ru 


Fe 


1500 


THF 


Br 


Pt 


Mn 


150 


DMFA 


Cl 


Ru 


Mn 


150 


DMFA 


Cl 


Pt 


Fe 


150 


DMFA 


Cl 


Ru 


Fe 


150 


DMFA 


Cl 


Pt 


Mn 


1500 


DMFA 


Cl 



1500 100 3360 

1500 100 2680 

1500 100 2640 

1500 100 3520 

1500 100 3620 

1500 100 3620 

1500 100 3830 

1500 100 2710 

1500 100 2580 

1500 100 3310 

1500 100 3740 

1500 100 2460 

1500 100 2410 

1500 100 2020 

1500 100 2300 

1500 100 3270 

1500 100 3250 

1500 100 2800 

1500 100 2520 

1500 100 3850 

1500 100 3570 

1500 100 3540 

1500 100 3950 

1500 100 2730 

1500 100 2520 

1000 120 3300 

1000 120 3800 

1000 120 2640 

1000 120 2520 

1000 120 2220 

1000 120 2270 

1000 120 3960 

1000 120 3360 

1000 120 2560 

1000 120 2440 

1000 120 3560 

1000 120 3620 

1000 120 3650 

1000 120 3570 

1000 120 2650 

1000 120 2710 

1000 120 3380 

1000 120 3430 

1000 120 2480 

1000 120 2310 

1000 120 2510 
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Ru Mn 1500 DMFA CI 1000 120 2410 

p t Fe 1500 DMFA CI 1000 120 3320 

» u Fe 1500 DMFA CI 1000 120 3620 

p t Mn 150 THF CI 1000 120 2880 

ru Mn 150 THF CI 1000 120 2430 

m pe 150 THF CI 1000 120 3940 

£ Fe 150 THF CI 1000 120 4290 

Pt Mn 1500 THF CI 1000 120 3730 

R u Mn 1500 THF CI 1000 120 3730 

p" pe 1500 THF CI 1000 120 2730 

ru Fe 1500 THF CI 1000 120 2480 

p7 Mn 150 DMFA Br 1500 120 3510 

R u Mn 150 DMFA Br 1500 120 3420 

p" pe 150 DMFA Br 1500 120 2340 

Ru pe 150 DMFA Br 1500 120 2390 

P? Mn 1500 DMFA Br 1500 120 2210 

Ru Mn 1500 DMFA Br 1500 120 2660 

Pt Fe 1500 DMFA Br 1500 120 3270 

ru Fe 1500 DMFA Br 1500 120 3550 

p7 Mn 150 THF Br 1500 120 2450 

ru Mn 150 THF Br 1500 120 2830 

J? Fe 150 THF Br 1500 120 3910 

ru pe 150 THF Br 1500 120 3850 

P? Mn 1500 THF Br 1500 120 3470 

ru Mn 1500 THF Br 1500 120 3730 

Fe 1500 THF Br 1500 120 2360 

ru Fe 1500 THF Br 1500 120 2660 

p t Mn 150 DMFA CI 1500 120 3640 

ru Mn 150 DMFA CI 1500 120 3440 

p" pe 150 DMFA CI 1500 120 2380 

Ru Fe 150 DMFA CI 1500 120 2230 

Pt Mn 1500 DMFA CI 1500 120 2510 

Ru Mn 1500 DMFA CI 1500 120 2150 

P? pe 1500 DMFA CI 1500 120 3100 

Z pe 1500 DMFA CI 1500 120 3810 

Pt Mn 150 THF CI 1500 120 2380 

ru Mn 150 THF CI 1500 120 2900 

KU THF CI 1500 120 3740 



Pt Fe 150 



ru Fe 150 THF CI 1500 120 3620 

p t Mn 1500 THF CI 1500 120 3610 

ru Mn 1500 THF CI 1500 120 3680 

Pt pe 1500 THF CI 1500 120 2570 

ru Fe 1500 THF CI 1500 120 2880 
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[000 1 ] Analysis of the data using matrix estimation formula (I) gives the 

information of TABLE 4. The Terms are the main effects and interactions of the 
factors in TABLE 3 and A-F are as given in the heading of TABLE 3. Thus A is the 
main effect of Factor "Primary Catalyst" and AD is the interaction effect of the 
Factors "Primary Catalyst" and "Cosolvent." Effects are p s. 
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Table 4 



Term 
A 
B 
C 
D 
£ 
F 
G 

AB 

AC 

AD 

AE 

AF 

AG 

BC 

BD 

BE 

BF 

BG 

CD 

CE 

CF 

CG 

DE 

DF 

DG 

EF 

EG 

FG 

ABC 

ABD 

ABE 

ABF 

ABG 

ACD 

ACE 

ACF 

ACG 

ADE 

ADF 

ADG 

AEF 

AEG 

AFG 

BCD 

BCE 

BCF 

BCG 

BDE 

BDF 

BDG 

BEF 



Effect 

-25 
2 
-43 
164 

32 

20 
229 

10 
2 
7 

22 
-26 
42 
39 
29 
73 
16 

-12 ^ 

6 
-28 

23 
-27 

28 
-31 
-35 

-9 
-25 
-188 
-21 

13 

-6 

4 
1 

-14 
1 

27 
-20 
30 
-17 
-6 
-53 
17 
53 
-1005 
-36 
14 
28 
-10 
47 
7 
36 



Term 

BFG 

CDE 

CDF 

CDG 

CEF 

CEG 

CFG 

DEF 

DEG 

DFG 

EFG 

ABCD 

ABCE 

ABCF 

ABCG 

ABDE 

ABDF 

ABDG 

ABEF 

ABEG 

ABFG 

ACDE 

ACDF 

ACDG 

ACEF 

ACEG 

ACFG 

ADEF 

ADEG 

ADFG 

AEFG 

BCDE 

BCDF 

BCDG 

BCEF 

BCEG 

BCFG 

BDEF 

BDEG 

BDFG 

BEFG 

CDEF 

CDEG 

CDFG 

CEFG 

DEFG 

ABCDE 

ABCDF 

ABCDG 

ABCEF 

ABCEG 



Effect 

22 
-24. 
-14. 
1 

-67 
20 
13 
1 

-16 
40 ; 

-373 
2 

-52 

-40- 
30 
44 v 
13 
21 

-U_ 
13 ■ 
20 
-4 • 

-22' 
12 . 

-40 
21 

-32 
7 

-41 
-1 
4 

-56 
14 

24 , 
-8 ■ 

29 : 

6 
6 

12 : 
-32 
30 
19 
18 
7 

58 , 
-2 - 
-8 - 
16 * 
29 

-23 

-34 



Term 

ABDEF 

ABDEG 

ABDFG 

ABEFG 

ACDEF 

ACDEG 

ACDFG 

ACEFG 

ADEFG 

BCDEF 

BCDEG 

BCDFG 

BCEFG 

BDEFG 

CDEFG 

ABCDEF 

ABCDEG 

ABCDFG 

ABCEFG 

ABDEFG 

ACD EFG 

BCDEFG 

AB CDEFG 



Effect 

-24 
-7 
-58 
9 
-25 
22 
-16 
51 
-6 
3 
21 
-40 
9 
-26 
18 
51 
1 
20 
-73 
-57 
-4 
4 
-11 



[000 1] The Effects are fitted to a normal probability plot and four points are 

identified as falling outside two standard deviations of the straight line fit: D, G, FG, 
and BCD. The BCD interaction is identified as a potential lead to nonlinear behavior 
simultaneously involving the Cocatalyst Metal (B), the Cocatalyst Metal Amount (C), 
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and the Cosolvent (D). Repeated followup iterations identify a strong synergistic 
effect of high levels of Fe in the presence of THF. 

[0001] While preferred embodiments of the invention have been described, the 

present invention is capable of variation and modification and therefore should not be 
limited to the precise details of the EXAMPLE. The invention includes changes and 
alterations that fall within the purview of the following claims. 
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